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Abstract 



We study the ground state and the collective excitations of parabolically- 
confined double-layer quantum dot systems in a strong magnetic field. We 
identify parameter regimes where electrons form maximum density droplet 
states, quantum-dot analogs of the incompressible states of the bulk integer 
quantum Hall effect. In these regimes the Hartree-Fock approximation and 
the time-dependent Hartree-Fock approximations can be used to describe the 
ground state and collective excitations respectively. We comment on the 
relationship between edge excitations of dots and edge magneto-plasmon ex- 
citations of bulk double-layer systems. 
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I. INTRODUCTION 



In recent years, a large body of work has been denoted to the study of double-layer 
two-dimensional electron systems! and also to the understanding of electronic properties 
of two-dimensional quantum dot systems. I In both cases such a huge interest has been 
triggered by advances in nanofabrication technology which made possible the synthesis of 
these artificial systems. In the case of double-layer systems,! some of the most interesting 
novel physics has been uncovered in the strong magnetic field limit of quantum Hall effect. 
Recent experimentafl and theoreticaliHi work in such strong magnetic fields has identified 
an unusual type of broken symmetry responsible for an unexpected quantum Hall effect 
which occurs at the Landau level filling factor of the individual layers v = vt/2 = 1/2. This 
broken symmetry state has spontaneous interlayer phase coherence, i.e. phase coherence 
even in the absence of tunneling between the quantum wells. This effect is produced by the 
Coulomb interaction between electrons in different layers. For quantum dot systems it is 
also true that most of the research in this context has been focused on electronic properties 
in a strong magnetic field.0 Much of the work on quantum dots in this regime is related 
to the existence of maximum- density- droplet (MDD) states,0'Ei which are the quantum dot 
analogs of the incompressible states responsible for the quantum Hall effect in bulk systems, 
and to the edge reconstructions^ which occur when these states become unstable. The 
recent experimental realizationS of layered quantum dot systems adds to the motivation for 
theoretical studies of these systems. In this article we discuss MDD states of double-layer 
quantum dot systems. Complementary numerical exact diagonalization study of double- 
layer quantum dot systems in a strong magnetic field have appeared recentlylll and some 
preliminary results0 from the present study have been reported earlier. 

The paper is organized as follows. In Section II the model used here for two-dimensional 
double-layer quantum dot systems is presented and discussed. In Section HI the stability 
limits for the z/^ = 1 MDD state are analyzed. This state is the quantum dot analog 
of the phase-coherent incompressible state in the bulk limit. In Section IV we discuss 
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collective excitations of this state with an emphasis on the interplay between the gapless 
edge excitations, alwayslll associated with the quantum Hall effect, and the gapless Goldstone 
modes of the broken-symmetry ground state in the bulk. In Section V We turn our attention 
to the ut = 2 MDD state , which corresponds to the bulk state with a filled Landau level 
in each layer. Our emphasis here is on the relationship between low energy excitations of 
the quantum dot system and coupled edge magneto-plasmon modes of bulk double-layer 
systems. Our results are briefly summarized in Section VI. 



II. THE MODEL 

We model a double-layer system made out of two-dimensional quantum dots by assuming 
identical parabolic potentials, V{r) = im*f2^r^, in the two layers. (In Section V it is 
emphasized that many of our results apply equally well to systems with bulk two-dimensional 
electrons in each layer.) The system is placed in a strong magnetic field B perpendicular to 
the layers. We restrict our attention to the strong magnetic field limit where Q/uJc << 1, 
and only the states in the lowest Landau level are important. (Here Uc = eB/m*c is the 
cyclotron frequency.) The lowest-Landau-level single-particle eigenstates in the symmetric 
gaugei are labeled by the angular momentum m: 

and 

Em = ^fiuJc + l{m + 1) (2) 

where 7 = m*f^^^^ = hud^/ujcy , i'^ = hc/eB, f = {x,y), z = x + iy is the 2D electron 
coordinate expressed as a complex number, and the allowed values of the single-particle 
angular momentum within the lowest Landau level are m = 0,1, 2,.... We will assume that 
the electron system is completely spin-polarized by the magnetic field. 

Including a phenomenological term describing tunneling between the two quantum dots 
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and up to an irrelevant constant, the second-quantized Hamiltonian of the system is given 
by 

"mu rnau' 

~'~2 ^ Kri'jmi,mim2'^m.'^(T'^m'2(T''"'Ti2o-'Cr7iio- (3) 

aa' 

where 



K^fm>im2 = {m[m'^\Vo\mim2) + aa' {m[m'^\V,\mim2) , (4) 

and t is the hopping amphtude between the two layers. The layer index is a =t, i where 
invoking a helpful§ analogy between the layer degree of freedom and the spin degree of 
freedom, cr =| corresponds to electrons in the right layer and a =i to electrons in the left 
layer. In Eq. H, Vq = {Va + Ve)/2 and = (Va — Ve)/2, are proportional to the sum and 
difference of intra-layer (Va) and interlayer (Ve) Coulomb interactions. 

For ut = 1 the ground state of the quantum dot in the Hartree-Fock approximation has 
its pseudospin polarized in the x — y pseudospin plane. Physically, electrons in this state 
occupy states which are a coherent linear combination of states localized in the separate 
layers. In order to study this state it is convenient to transform to a new representation, 
carrying out a rotation in pseudospin space by defining 

cL = -^{(^(^La + "L)- (5) 

In this representation, a pseudospin up electron is in a symmetric double-layer state and a 
pseudospin down electron is in an antisymmetric state. The Hamiltonian can alternately be 
expressed in the form: 

n = ^(m7 - at)al^„ama + 

ma 

r?iim2m,3m4(T(T' 

where a =|, | for a =t, |. In this rotated pseudospin representation the hopping parameter 
t simply acts as an external magnetic field. We will discuss here only the limit t — > 0. 



III. THE STABILITY OF MDD STATE 



For non-interacting electrons in a single-layer quantum dot, the many-body ground 
state is a single Slater determinant in which the confinement energy is minimized by occupy- 
ing the orbitals from m = Otom = A^— 1. This state is an exact many-body eigenstate of the 
Hamiltonian even when electron-electron interactions are included.0 We refer to this state as 
the maximum density droplet (MDD) state. The MDD state is the analog for quantum dots 
of the bulk u = 1 state. For double layers a new type of MDD state can occur which is the 
analog of the bulk double-layer ut = 1 broken-symmetrylil ground state.ii'iil In this state 
electrons occupy symmetric states (pseudospin up) with m = Otom = A^ — 1 such that the 
total angular momentum is Mq = N{N — l)/2. Many-body states with smaller total angular 
momentum than this MDD state will have the advantage of smaller confinement energies 
since the electrons are closer to the minimum of the confinement potential, but the disad- 
vantage of larger interaction energies since the electrons are closer to each other. For weak 
confinement the angular momentum of the ground state of an interacting-electron droplet 
will occur at M > Mq in order to reduce the Coulomb energy while for sufficiently strong 
confinement the ground state will occur at M < Mq in order to reduce the confinement 
energy. 

A necessary condition for the MDD state to be the ground state is that all occupied or- 
bitals have a Hartree-Fock quasiparticle energy which is lower than the quasiparticle energies 
of all unoccupied orbitals. Using Eq. (§) we have found that the Hartree-Fock quasiparticle 
energies for the symmetric MDD state are given by 

to' to' 

where 

^mm' — {fTim\Vo\mm') — {m'm\Vo\mm'), (8) 

^mm' — {'fnm'\Vo\mm') — {m! m\Vz\mm') . (9) 
For the symmetric MDD state Uma is 1 for cr =1 and m < N, and it is zero otherwise. 
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FIG. 1. (a) The single particle H-F energy for the double-layer quantum system with 20 elec- 
trons and interlayer distance d = £ at different confinement potential strengths, (b) The phase 
diagram for MDD states in = 20 double-layer quantum dots. The ground state in the region 
of "compressible" state tends to have many low energy excitation away from the edge of quantum 
dots. 



Fig. |T]a shows Hartree-Fock quasiparticle energies for symmetric and antisymmetric or- 
bitals for = 20, t = 0, and d = i at three different confinement strengths. The exchange 
splitting between the occupied symmetric orbitals and the unoccupied antisymmetric or- 
bitals is due to interlayer interactions and decreases in magnitude as d increases or as m 
increases toward the edge of the dot. For 7/(e^/e£) = 0.06 the confinement potential is 
relatively weak and unoccupied symmetric orbitals near the edge of the dot are lower in 
energy than occupied symmetric orbitals near the center of the dot. For 7/(e^/e£) = 0.12 
the confinement potential is now strong and unoccupied antisymmetric orbitals near the 
center of the dot are lower in energy than occupied symmetric orbitals near the edge of the 
dot. These results show that there is a finite interval of confinement strengths, including the 
case 7/(e^/e£) = 0.10, for which the MDD state is stable. Fig. [pj shows a phase diagram 
constructed for = 20 using this quasiparticle energy stability criterion, which shows how 
the interval where the confinement potential strength renders the MDD state stable nar- 



6 



rows and shifts as a function of layer separation. (The actual stability region for the MDD 
state beyond the Hartree-Fock approximation will be slightly narrower because of excitonic 
corrections to particle-hole excitation energies). For d —>■ interlayer and intralayer interac- 
tions become identical and these results would apply equally well to the ut = 1 state for a 
single-layer system in the limit of vanishing Zeeman coupling (Note that we assume strong 
Zeeman coupling for the double-layer systems discussed in this paper). For larger the 
region of stability of a parabolically confined MDD droplet becomes smaller; for very large 
the electron density in a parabolic external potential must follow the semi-elliptic behav- 
ior preferedS by electrostatics rather than the nearly constant charge density profile of the 
MDD state. For sufficiently large N the MDD state is never stable in a parabolic external 
potential and can be stabilized only by the potential from an approximately neutralizing 
positively charged background. 



IV. COLLECTIVE EXCITATIONS 

In the time-dependent Hartree-Fock approximation, elementary excitations of the MDD 
state are constructed by allowing its particle-hole excitations to couple. The Hamiltonian for 
the double-layer system does not mix states in which the number of particles in symmetric 
(pseudospin |) states differs by an odd integer; elementary excitations with even and odd 
numbers of pseudospin reversed states do not mix. In order to describe the odd excitations 
we introduce the operators 



PM,in ~ '^■m+M,l 



am,] and ptM,m — Oim,](^m-M,i- 



(10) 



The equations of motion for these operators are readily derived and reduce to a closed system 
of equations after Hartree-Fock factorization: 



[n,p'. 



■M,m. 



J 



mm' mm' 



I 



\ 



-E. 



-M 



Pm. 



(11) 



mm' j y P—M,m' j 



where 



{{m,m' + M\V^\m + M,m') - {m' + M,m\Vo\m + M,m')) (12) 



and 



-T^^; = {mm \V^\m + M,m' - M) - {m' m\Vz\m + M , m' — M) . (13) 

Diagonalizing the above matrix will give us the time- dependent Hartree-Fock approxima- 
tion (TDHFA) pseudospin-flip elementary excitations of the MDD state. If these operators 
evolved according to the mean-field Hartree-Fock Hamiltonian for the ground state, the 
excitation energies would simply equal the difference of occupied and unoccupied Hartree- 



Fock eigenvalues. The additional terms in Eq. ([Tl| ) reflect the changes in the mean-field 
Hamiltonian in the excited state and, in a diagrammatic derivation of the TDHFA , would 
arise from vertex and repeated-bubble corrections to the single-loop approximation for the 
pseudospin-flip response function. Below we refer to the influence of these additional terms 
as vertex corrections. 

The even elementary excitations are found by considering the equation of motion for the 
operator p% jn = Q^m+M.t'-'^ni.T- Hartree-Fock factorization of the equation of motion of these 
operators again gives a closed set of equations: 

Diagonalizing the matrix G^^, gives the even elementary excitations of the MDD state. In 
Eq. (0) 

(^M _ r / HF _ f.HF\ 
^mm' ~ "mm' l^^m+Aft mj / 

+ {{mm' + M\Vo\m + Mm') - {m' + Mm\Vo\m + Mm')). (15) 

Again vertex corrections cause the excitation energies to differ from the difference of Hartree- 
Fock eigenvalues. 
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M 

FIG. 2. Isospin-flip excitations at vt = ^ with vertex correction and without vertex correction 
for a double-layer quantum dot system with = 20 electrons, interlayer distance d = £, and 
confinement potential 7 = O.lOe^ /£. 

Fig. ^ shows the odd elementary excitations obtained for layer separation d = i, N = 20, 
and confinement potential 7 = 0.10(e^/£) as a function of angular momentum M. Note that 
for each M the excitation energy has a contribution M7 from the confinement potential. The 
number of particle-hole excitations with angular momentum change M is for M > and 
N+M for M < 0. (The minority-spin angular momentum cannot be negative.) When vertex 
corrections are neglected, the lowest excitation energies occur at the most negative angular 
momentum and correspond to transferring an electron from pseudospin up at the edge of 
the quantum dot to pseudospin down at the center of the quantum dot. Vertex corrections 
change the elementary excitation spectrum qualitatively. The most obvious modification 



9 



is the appearance of a zero-energy excitation for M = which we discuss below. The 
M = —N + 1 excitation energy is also reduced by vertex corrections; with increasing 7 
the energy of this excitation will become negative before the quasiparticle energy stability 
criterion is violated so that the region of stability of the MDD state is narrower than implied 
by this criterion. A similar reduction in the excitation energies occurs at positive M. 

In bulk systems the broken symmetry ground state of the ut = 1 double-layer system 
has a gapless Goldstone modei with linear dispersion (It can be shown that the M = zero- 
energy elementary excitation just corresponds to a global rotation of the pseudospin). Some 
remnant of these Goldstone modes should be present in the excitation spectrum of finite-size 
quantum dot systems. In the limit of large dots, we expect that M/ R should act like the 
azimuthal component of a two-dimensional wave vector and that the many collective modes 
that occur at a given M are approximately related to the discrete set of radial wavevectors 
geometrically defined by the finite radius Rn ~ V2N£ of the quantum dot which is related, 
at least approximately, to the zeroes of Jo^qR^)- We have, however, not yet been able to find 
a completely satisfactory semiclassical interpretation of our microscopic TDHFA results in 
this way. For example the pronounced asymmetry between the low-energy TDHFA collective 
modes for M < and M > 0. 

The even elementary excitations represent charge-density-wave edge excitations of the 
incompressible MDD state and occur only for M > 0. In Fig. ^ we show TDHFA re- 
sults calculated for = 20, 7 = 0.1(e^/e£), and d = £. Again the excitation energies 
are substantially reduced by vertex corrections. These results are identical to those that 
would be obtained for a single-layer system with an interaction which is the average of the 
double-layer intra-layer and inter-layer interactions. As discussed in detail previously for the 
single-layer case,0 the M = 1 excitation energy is exactly equal to 7 and is not influenced 
by electron-electron interactions. These charge density wave excitations of the quantum dot 
are harbingers of the edge magneto-plasmon excitations of bulk systems which we discuss 
at greater length in the following section. With decreasing confinement strength a collective 
mode will become unstable slightly before the quasiparticle-energy stability criterion is vio- 
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lated, narrowing the region of stability of the MDD state on the weak confinement side just 
as it is narrowed on the strong confinement side. Very recently^ for the d = case, Sondhi 
et al. have identified topological edge instabilities which further limit the stability region of 
MDD states. 
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FIG. 3. Edge Charge density-wave edge excitations with vertex correction and without vertex 
correction at i/j- = 1 for a double-layer quantum dot system with = 20 electrons and interlayer 
distance d = £ and confinement potential 7 = 0.10(e^/^). 



V. EDGE MAGNETO-PLASMONS AND z^t = 2 DOUBLE-LAYER QUANTUM 

DOTS 

The quantum dot analog of the bulk v = 2 state in a double-layer system is the single 
Slater determinant state in which the single-particle angular momentum states from m = 
to m = — 1 are occupied in each layer. It can be shown that this state is a pseudospin 
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singlet. In quantum dot systems this state becomes unstable with increasing confinement 
strength only when it becomes energetically favorable to occupy higher Landau levels, and 
becomes unstable with decreasing confinement strengths when it becomes favorable to move 
electrons away from the center of the quantum dot, for example by edge reconstruction. 
For parabolic confinement, the region of stability depends on d and narrows with increasing 
N. The excitation spectrum of this state is less subtle than for the v = 1 case discussed 
above because a gap ~ fiujc exists in the bulk for both pseudospin-fiip and charge density 
excitations. The only low-energy excitations are locahzed at the edge of the system. 

In this section we work in the representation where pseudospin eigenstates are localized 
in individual layers. We consider the equations of motion for the operators which create in 
phase and out-of-phase density wave excitations at the edges of the two-layers: 

Pm,M — ^m+M,a^m,a (16) 

and 

Pm,M — '^'^(^m+M,a^m,a, (17) 

where m < N and m + M > N . After a Hartree-Fock factorization the equations of motion 
for these operators close and a calculation very similar to that detailed in Section IV gives 
the following result: 

["^'Pm.M] — {^m+M ~ ^m^)Pm,M (18) 

+ I3Pm',M[2(^,^' + M|Vo|m + M,m') - (m, m' + M| V^jm', m + M)] (19) 

m' 

and 

["^^ Pm,M] — {^m+M ~ ^m^)Pm,M (20) 

+ I]Pm',M[2(^,^' + M|l^^|m + M,m') - (m, m' + M| Vk|m', m + M)]. (21) 

m' 

In this case the single-particle Hartree-Fock energies, 

— 7(m -|- 1) -|- ^ [2(^71, m'|Vo|m, m') — (m, m'| VOi|m', m)], (22) 

m' 
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are pseudospin independent. 



2.0 




0.5 




0.0 

0.0 5.0 10.0 15.0 20.0 

M 

FIG. 4. Spin-density wave edge excitations at z^j^ = 2 for a double-layer system with = 20, 
interlayer distance d = i, and confinement potential 7 = 0.20(e^/^). 

Illustrative TDHFA quantum dot collective excitation energies at ut = 2, evaluated using 
these expressions, are shown in Figs. ^ and ^ We comment primarily on the M = 1 results 
which we relate to edge magneto-plasmons of bulk double-layer systems below. For M = 1, 
m = A^ — 1 in both Eq. (|1^) and Eq. (|T7|), explicit expressions can be given for the excitation 
energies: 

N-1 

Ef^=7+ J2 [2{N,m'\Vo\N,m) - {N,m'\VA\m',N)] 

m'=0 

N-1 

- J2 [2(^^-1, m'\Vo\N - 1, m') - {N - 1, m'\VA\m', N - 1)] 

m'=0 

+2{N - 1, N\Vo\N, N - 1) - {N - 1, N\Va\N - 1, N), (23) 
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and 

N-1 

e{^)=7+ J2 [2{N,m'\Vo\N,m') - {N,m'\VA\m',N)] 

m'=0 
N-1 

- J2 [2(A^ - 1, m'\Vo\N - 1, m) -{N -I, m'\VA\m\ N - 1)] 

m'=0 

+2{N - 1, iV|V;|A^, - 1) - (A^ - 1, N\Va\N - 1, A^). (24) 
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FIG. 5. Charge-density wave edge excitations = 2 for a double-layer system with = 20, 

interlayer distance d = £, and confinement potential 7 = 0.20(e^/^). 

The in-phase mode has interaction corrections from the difference of the interaction 
contribution to the Hartree-Fock single-particle energies at m = A^ and m = A^ — 1 and 
also from vertex corrections. In the case of single-layer quantum dots0 these two classes of 
corrections are known to cancel exactly. The easiest way to establish this property in the 
single layer case is to note that the M = 1 excitation operator simply changes the state 
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of the center-of-mass of all electrons without changing the relative motion state to which 
interactions are sensitive. For Ei'^ the same argument goes through in the double layer case 
and only the external potential contribution to the excitation energy survives: 

= 7- (25) 

This cancelation is verified numerically as seen in Fig. g Eq. (H) can also be derived from 
the following identities which are established for interaction matrix elements in the lowest 
Landau level: 

Af-l N-l 

{N,m'\V\N,m') - ^ (A^ - 1, m'|V"|Ar - 1, m') = -(A^ - 1, A^|1/|iV, iV - 1) (26) 

m'=0 m'=0 

and 

N-l N-l 

Y {N,m'\V\m',N) - ^ (A^ - 1, m'|1/|m', iV - 1) = -(A^ - 1, iV|V"|Ar - 1, iV). (27) 

m'=0 m'=0 

These identities hold for any interaction potential and in particular for = Vq or 1^ = V^. 
It follows from Eq. (g§) and Eq. d^^) that 

e['^ =-f -2{N -1,N\Ve\N,N -1). (28) 

This simple result has also been verified numerically. 

These M = 1 edge excitations are the quantum dot precursors of the edge magneto- 
plasmon collective excitationsH of larger A^ electronic systems in strong magnetic fields 
which give rise to far-infrared resonances in micron scale systems and to surprisingly sharp 
radio-frequency resonances in large area 2DEG's. Some helpful remarks on edge magneto- 
plasmons in double-layer systems follow from the above analysis. As we have emphasized, 
larger A^ MDD states are bulk quantum Hall states and they will occur only when the ex- 
ternal potential is quite close to that from a neutralizing background of positive charge. For 
a background charge density which neutralizes the electron density in each layer it follows 
from Eq. (^) that the difference between the expectation values of the external potential 
for the states with angular momenta A^ and A^ — 1 is 
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e2 



j = {n-i,n\Va + Ve\n, slHR/^) + HR/d)]. (29) 

The final form for 7 assumes an electron disk whose radius R is large compared to both 
the magnetic length i and the layer separation d. (We assume that d is at least ~ £.) 
The two terms in square brackets in Eq. (EUI) comes respectively from intra-layer and inter- 
layer interactions with the background. As discussed above, the energy of the in-phase 
collective magnetoplasmon excitations of the double-layer system (-E^pi) depends only on 
the external potential and is therefore given by E^^i = 7. Since (A^ — 1, Ve|A^, — 1) ~ 
(e^/vri?) \ia{R/d) for d it follows that the energy of the out-of-phase mode is 

E-^, = ^[HR/i)-ln{R/d)]. (30) 

For large 2DEG's the splitting of the edge magnetoplasmon energies due to coupling between 
the layers persists to very large layer separations; for example E^^^i = 'iE^^i for d = (RiY^"^. 
When the separation between the two layers is comparable to i, often a requirement for 
observable coupling effects in bulk systems, the edge magnetoplasmon coupling will typically 
make E^^i unobservably small. We remark that these results have been derived for disorder 
free systems. We expect that these results for edge magnetoplasmon energies will remain 
valid in disordered systems provided that the microscopic length scale i, which describes the 
degree of localization of the edge wave in the clean limit, is replaced by the appropriatelHl 
disorder-dependent length. 



VI. SUMMARY 

We have studied double-layer quantum dot systems focusing on the MDD droplet states 
corresponding to the bulk incompressible quantum Hall states at z/ = 1 and u = 2. Ground 
states and collective excitations have been approximated using Hartree-Fock and time- 
dependent-Hartree-Fock approximations. The Hartree-Fock u = 1 MDD ground states 
have only symmetric double-layer orbitals occupied and have spontaneous interlayer phase 
coherence. The regime of stability of this state has been estimated as a function of layer 
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separation and confinement potential strength. The broken symmetry in this state gives 
rise to low-energy excitations geometrically confined 'Goldstone' collective excitations of the 
quantum dot. The Hartree-Fock u = 2 MDD ground state has orbitals in both layers occu- 
pied. The only low-energy collective excitations of this system are locahzed near the edge 
of the system and correspond to coupled magnetoplasmon excitations of bulk double-layer 
systems. We have used our microscopic calculations for quantum dots to derive expressions 
for the magnetoplasmon energies of bulk systems. 
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